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Zero- and longitudinal-field muon spin relaxation (/iSR) experiments have been carried out in the 
alloy series Pr(Osi_ I Ru 2! )4Sbi2 and Pri_ y La H Os4Sbi2 to elucidate the anomalous dynamic muon 
spin relaxation observed in these materials. The damping rate A associated with this relaxation 
varies with temperature, applied magnetic field, and dopant concentrations x and y in a manner 
consistent with the "hyperfine enhancement" of 141 Pr nuclear spins first discussed by Bleaney in 
1973. This mechanism arises from Van Vleck-like admixture of magnetic Pr 3+ crystalline-electric- 
field-split excited states into the nonmagnetic singlet ground state by the nuclear hyperfine coupling, 
thereby increasing the strengths of spin-spin interactions between 141 Pr and muon spins and within 
the 141 Pr spin system. We find qualitative agreement with this scenario, and conclude that electronic 
spin fluctuations are not directly involved in the dynamic muon spin relaxation. 

PACS numbers: 71.27. +a, 74.70.Tx, 74.25.Nf, 75.30.Mb, 76.75.+i 



I. INTRODUCTION 

The filled skutterudite compound PrOs4Sbi2 is the 
first praseodymium-based heavy-fermion superconduc- 
tor to be discovered, 1 and is one of the few /-electron 
heavy-fermion compounds in which a non-Kramers ion 
exhibits a nonmagnetic crystalline-electric-field (CEF)- 
split ground state. Both the normal and superconduct- 
ing states of PrOs4Sbi2 are unusual: The T\ singlet 
Pr 3 + ground state is separated from a first excited 
state (tetrahedral notation 2 ) by a remarkably small split- 
ting AScefAb « 7 K, 3 ~ 6 leading to strong CEF ef- 
fects at low temperatures. There is no sign of mag- 



netic ordering. The Sommerfeld coefficient 7 is diffi- 
cult to measure in the presence of the CEF Schottky 
anomaly in the low-temperature specific heat, but is 
estimated to lie between 500 and 750 mJ mol -1 K -2 . 
The Pr 3+ ions are enclosed in an icosahcdral cage of 
Sb atoms that is considerably larger than the ionic 
size, and large-amplitude Einstein-like 4/-ion phonon 
modes ("rattling" modes) are observed. Below the su- 
perconducting transition temperature T c = 1.85 K an 
unconventional superconducting phase is found, 7 with 
evidence for multiple phases, time reversal symmetry 
breaking, 8 and extreme multiband behavior. 9 ' 10 Disper- 
sive antiferroquadrupolar excitons and an unusual low- 
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temperature high-field phase with antifcrroquadrupolar 
order are observed. The mechanism or mechanisms for 
heavy-fermion behavior and Cooper pairing in this com- 
pound remain controversial; rattling modes 11 and Pr 3+ 
quadrupole fluctuations 12 have been proposed. 

Muon spin relaxation (^iSR), 13,14 like other magnetic 
resonance techniques, probes magnetism and electronic 
structure in solids on the microscopic (atomic) size scale. 
In fiSK experiments spin-polarized positive muons 
are implanted into the material of interest and stop at 
interstitial sites. During its lifetime each muon spin pre- 
cesses in the magnetic field at its site and eventually de- 
cays — > e + + v e + V^); the direction of the emitted 
positron's momentum is correlated with the muon spin 
direction at the moment of decay. A large number of such 
events determines the time development (relaxation) of 
the ensemble-average muon spin polarization. Several 
fiSR studies of Pr-based filled skutterudites have been 
reported. 8 ' 15 ^ 20 

The magnetic environment of the muon creates a local 
field Hioc at the muon site that causes muon spin relax- 
ation. Relaxation mechanisms can be divided into two 
classes, depending on the behavior of Hi oc : 

• static (or quasistatic 21 ) relaxation, due to an inho- 
mogeneous static distribution of H loc that causes a 
spread of muon Larmor frequencies and consequent 
loss of muon spin phase coherence; 

• dynamic relaxation (often called spin-lattice relax- 
ation in the NMR literature), due to thermal fluc- 
tuations of Hioc that induce transitions between 
muon spin levels and equilibrate the muon spin 
populations. The equilibrium muon spin polariza- 
tion is negligible compared to the initial polariza- 
tion (~100%). The dynamic component of Hi oc 
usually arises from electronic spin fluctuations; we 
shall see, however, that nuclear spin fluctuations 
can also be involved. 

In longitudinal- field ,uSR (LF-^iSR) , which includes zero- 
field fiSR (ZF-/xSR) as a special case, a magnetic field Hl 
is applied parallel to the initial muon spin polarization. 
The dependence of the muon spin relaxation on Hl helps 
to separate the static and dynamic contributions to the 
relaxation. 22,23 

The ZF-/iSR spin relaxation function in PrOs4Sbi2 at 
low temperatures 8 could be fit by the product of an expo- 
nential damping factor cxp(— At) and the Kubo-Toyabe 
(K-T) functional form 22 ' 23 expected for a Gaussian qua- 
sistatic field distribution. This behavior was attributed 
to a two-component form of Hi oc : a static component, 
responsible for the K-T relaxation, and a second compo- 
nent, responsible for the exponential damping, which was 
determined by LF-^SR measurements to be due to dy- 
namic fluctuations. 8 The origin of this dynamic compo- 
nent was not clear, although 4/-clcctron dynamics asso- 
ciated with the small CEF splitting were noted as a pos- 
sible mechanism. In the normal state above T c the qua- 



sistatic field distribution was attributed to dipolar inter- 
actions between muon spins and neighboring (principally 
121 Sb and 123 Sb) nuclear magnetic moments. An increase 
in the quasistatic relaxation rate was observed below T c , 
and interpreted as evidence for time-reversal symmetry 
breaking in the superconducting state of PrOs4Sbi2. 8 
This effect has also been studied in Pr(Osi_a;Ru x )4Sbi2 
and Pri_ y La y Os4Sbi2 alloys, 24,25 which are supercon- 
ducting for all x and y with T c of the order of 1 K. 26,27 

In this article we report results of ZF- and 
LF-ytiSR experiments in the Pr(Osi_a;Ru x )4Sbi2 and 
Pr 1 _j / La J) Os4Sbi2 alloy systems, 28 which have been un- 
dertaken to elucidate the anomalous exponential damp- 
ing. We confirm that the damping is indeed dynamic 
in nature, and argue that it is due to enhancement 
of 141 Pr nuclear magnetism via intra-atomic hyperfinc 
coupling to the Pr 3+ 4/ electrons. 29 This coupling in- 
duces a Van Vleck-like admixture of low-lying Pr 3+ mag- 
netic CEF-split excited states into the singlet ground 
state, thereby increasing the effective 141 Pr nuclear mo- 
ment. Both the 141 Pr-/i + and the 141 Pr- 141 Pr dipolar 
interactions are increased (by factors of ^20 and ^400, 
respectively) by this "hyperfinc enhancement" mecha- 
nism, which also increases indirect 141 p r _ 141 p r inter- 
actions mediated by exchange coupling between Pr 3+ 
ions. 29 Hyperfine-enhanced 141 Pr nuclear moments have 
a significant effect on muon spin relaxation: the mecha- 
nism accounts for the observed damping of the muon spin 
relaxation, and is qualitatively consistent with the behav- 
ior of the damping rate A with enhancement strength and 
Pr-ion concentration dependence (for Pri_ 2/ La i) Os4Sbi2 
alloys) across both alloy series. Exchange-mediated in- 
teractions appear to dominate the 141 Pr spin dynam- 
ics. The origin of the dynamical muon spin relaxation 
has been clarified, and we conclude that electronic spin 
fluctuations (except those associated with the hyperfine- 
enhancement mechanism) are not directly involved in the 
anomalous muon spin relaxation. 

The article is organized as follows. After a brief 
description of the experiments in Sec. II, we report 
our experimental results for P^Osi-zRu^Sb^ and 
Pri_j / La y Os4Sbi2 in Sec. III. We consider the depen- 
dence of the muon exponential damping rate A on Hl, 
temperature, and dopant concentrations x and y, and 
their implications for the mechanism for the dynamic 
muon spin relaxation in Sec. IV. Our conclusions are 
summarized in Sec. V. 



II. EXPERIMENTS 

ZF- and LF-^iSR experiments were carried out in 
Pr(Osi- a; Ru ;E )4Sbi2, x = 0.05, 0.1, 0.2, 0.6, and 1.0, and 
(Pri_ y Laj / )Os 4 Sbi2, y = 0.2, 0.4, 0.6, and 0.8. ZF- and 
LF-/iSR experiments were carried out on powdered sam- 
ples at the Meson Science Laboratory, KEK, Tsukuba, 
Japan, and the MuSR spectrometer at the ISIS neu- 
tron and muon facility, Rutherford Appleton Laboratory, 
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Chilton, U.K. LF-/1SR experiments were carried out at 
the M15 beam line at TRIUMF, Vancouver, Canada, on 
mosaics of oriented ~1 mm 3 crystals, prepared by the 
Sb-flux method, 30 with (100) directions parallel to the 
applied field. 

3 He- 4 He dilution cryostats were used to obtain low 
temperatures. The mosaic crystals were mounted on a 
thin GaAs backing, which at low temperatures rapidly 
depolarizes muons and minimizes any spurious signal 
from muons that do not stop in the sample. Standard 
time-differential /iSR asymmetry data 13 ' 14 were taken in 
the normal and superconducting states at temperatures 
in the neighborhood of T c for Hl between and 125 Oe. 



III. EXPERIMENTAL RESULTS 

A. Zero-field muon spin relaxation 

Figure 1 shows the ZF-/iSR positron count rate 
asymmetry 13 ' 14 AG(t), where A is the initial asymmetry 
and G(t) is the muon spin polarization (initially ~100%), 
in Pr(Oso.gRuo.i)4Sbi2. The time-independent offset in 

0.2 
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FIG. 1: (color online) Zero-field positron count rate asymme- 
try (muon spin relaxation function) in Pr(Oso.9Ruo.i)4Sbi2, 
T = 1.31 K. Solid (blue) curve: fit of the exponentially 
damped Gaussian K-T function [Eq. (2)] to the data. Dashed 
(green) curve: best fit of a simple exponential to the data. 
Dash-dot line: time-independent signal from muons that 
stopped outside the sample. 



Fig. 1 comes from nonrelaxing muons that stopped in the 
silver cold finger or elsewhere in the cryostat (i.e., did not 
stop in either the sample or the GaAs backing). 

The zero-field static Gaussian K-T relaxation 
function 22,23 



f 



G* 1 (A, t) = - + |(1 - A 2 * 2 ) exp(-i A 2 * 2 ) (1) 

describes muon spin relaxation by randomly-oriented 
quasistatic muon local fields, with Cartesian components 
that vary according to a Gaussian distribution with a 



zero mean value and an rms width A / 7^ (7^ is the muon 
gyromagnetic ratio) . The "damped Gaussian K-T" func- 
tion 



G(t) = exp(-At) Gf l (A,t) 



(2) 



where A is the damping rate, was fit to the experimental 
data (solid curve in Fig. 1). Although the data appear 
to decay roughly exponentially, it can be seen that the 
best fit to a simple exponential (dashed curve in Fig. 1) 
is not as good as the fit to Eq. (2). 

Before analyzing the relaxation data further, we dis- 
cuss the choice of Eq. (2) as a fitting function. In ZF- 
and LF-/iSR, fluctuations of Hi oc arc often treated us- 
ing the "dynamic K-T" relaxation function, 23 in which 
Hi oc fluctuates as a whole with a single correlation time. 
Such a procedure is not appropriate in the present ex- 
periments, however. A quasistatic contribution to Hi oc , 
due mainly to Sb nuclear dipolar fields, is also present, 
so that Hioc is the sum of quasistatic and fluctuating 
components. 8 

Figure 2 shows ZF-/iSR data obtained from the 
isostructural compound LaOs4Sbi2 (here the time- 
independent signal has been subtracted), in which there 
are no 4/ electrons. 16 The best fit of Eq. (2) to these 
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FIG. 2: (color online) Zero-field muon spin polarization decay 
function in LaOs 4 Sbi 2 , T = 0.91 K. Data from Ref. 16. The 
signal from muons that do not stop in the sample has been 
subtracted. There is no significant exponential damping. 



data yields a negligible exponential damping rate, in- 
dicating that the damping arises from the presence of 
Pr ions. The values of the K-T relaxation rate in the 
normal state differ by only ~15% between PrOs4Sbi2 
[A = 0.143(5) lis' 1 } and LaOs 4 Sbi 2 [A = 0.167(5) ^s" 1 ]. 
In the latter compound the nuclear dipolar field is the 
only mechanism for the static K-T term. 16,31 

This near equality is additional evidence that A is due 
to nuclear dipolar fields in PrOs4Sbi2- Furthermore, sig- 
nificantly better fits to zero-field data are obtained using 
the damped static K-T function (reduced \ 2 typically 
~1.1) than the dynamic K-T function (reduced % 2 typ- 
ically ~1.4). We therefore use damped static relaxation 
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of the form of Eq. (2) to model the situation where muon 
spin states are split into Zeeman levels by the quasistatic 
component of + Hi oc , and the fluctuating component 
of Hi oc induces transitions between these Zccman levels. 

In Fig. 3 the damping rate A in zero field for 
Pr(Osi_ x Ru x ) 4 Sbi2, x = 0.1, 0.2, and 1.0, is plotted vs 
temperature T for 0.02 K < T < 2 K. As in PrOs 4 Sbi 2 , 8 



heat of PrOs4Sbi2, 3 is similarly enhanced. Hyperfinc- 
enhanced relaxation is independent of temperature for 
T <C AEcef/Izb, 34 which is roughly satisfied in the 
present experiments. 



B. Longitudinal-field muon spin relaxation 
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FIG. 3: Zero-field exponential damping rate A vs temperature 
in Pr(Osi_ a; Ru a; )4Sbi2, x = 0.1, 0.2, and 1.0. The vertical 
lines indicate the superconducting transition temperatures for 
each alloy. 



little temperature dependence is observed in this temper- 
ature range; this is also the case for Pri_ 2/ La J) Os4Sbi2 
alloys (data not shown). This weak temperature depen- 
dence strongly suggests that A is due to nuclear mag- 
netism rather than electronic spin fluctuations: the latter 
would be expected to show significant temperature de- 
pendence, especially below T c , whereas fluctuations aris- 
ing from nuclear spin-spin interactions are temperature- 
independent except at very low temperatures. 32 The dy- 
namic relaxation is unlikely to arise from fluctuations of 
121 Sb and 123 Sb nuclear spins, which should result in qua- 
sistatic contributions to Hi oc because their spin-spin and 
spin-lattice relaxation times are relatively long. 33 

A candidate mechanism for dynamic muon spin relax- 
ation is the hyperfine-enhanced 141 Pr nuclear spin system 
with hyperfine-enhanced effective nuclear moments, 29 , 
described briefly in Sec. I, which as discussed below in 
Sec. IIIB leads to dynamic muon spin relaxation rates 
in qualitative agreement with experiment. 141 p r nu- 
clear spin fluctuations due to hyperfine-enhanced 141 Pr- 
141 p r gpi n _ S pi n interactions can be rapid enough to cause 
dynamic muon spin relaxation, since both dipolar and 
exchange-mediated contributions to these interactions 
are increased by hypcrfine enhancement. The observed 
nuclear Schottky anomaly in the low-temperature specific 



We first review the effect of longitudinal applied 
field Hl on LF-/iSR measurements. 22,23 In the case where 
Hi oc is quasistatic, for Hl 3> A/7 M the (quasistatic) re- 
sultant field Hioc + Hi is nearly parallel to the initial 
muon spin direction. Then the precession that causes 
quasistatic relaxation is reduced in amplitude, and the 
muon spin polarization is nearly time-independent. For 
weaker fields Hl ~ A/7 M , the muon polarization at long 
times increases with increasing Hl- This phenomenon, 
called "decoupling," does not occur for dynamic relax- 
ation unless the latter is quenched for Hl ~ A/7 M ; this 
is unusual because nuclear dipolar fields are small, of the 
order of a few Oe, and applied fields this small seldom 
affect dynamic relaxation mechanisms. Thus LF-/1SR 
measurements help to determine whether the observed 
ZF relaxation is due to static or dynamic contributions to 
Hi oc . 23 In favorable cases the dependence of the damping 
rate A on Hl yields statistical properties (rms amplitude, 
correlation time) of the fluctuating field. 23 ' 35 

LF-^SR experiments were performed in weak longitu- 
dinal fields Hl in the normal state just above T c . Repre- 
sentative LF-/zSR spectra are shown in Fig. 4. The data 
exhibit the late-time field dependence that is a character- 
istic feature of decoupling, together with overall damping 
that is stronger at lower fields. A "damped static longi- 
tudinal K-T function" appropriate to nonzero Hl, of the 
form of Eq. (2) with G* T (A,i) now the static Gaussian 
K-T function in nonzero applied longitudinal field, 23 was 
therefore fit to the data. "Global" fits to all the field- 
dependent data at a given temperature were carried out, 
with A taken to be independent of field but varied for 
best fit. The field dependence of A was obtained under 
this condition. 

Figure 5 shows these field dependencies for several al- 
loys. The fit values of A are given in Table I, which 
also shows parameters derived from the field dependence 
of A as discussed below in Sec. IV A. The values of A 
obtained in this way are larger than those obtained from 
fits to ZF-/iSR data, and moreover A exhibits an increase 
with increasing Hl (not shown) if allowed to vary with 
field. This behavior indicates that the data depart some- 
what from the form of Eq. (2), and we consider the fit 
values of A and A to be of qualitative significance only. 

If the exponential damping were due to a (Lorentzian) 
distribution of static contributions to Hi oc , the observed 
ZF value of A leads to an estimate A/7 M ~ 1 Oe for the 
spread of the local fields H\ oc . Then an applied longi- 
tudinal field of order 10 Oe should nearly decouple H\ oc 
and there should be almost no damping. But it can be 
seen from Fig. 5 that in general A is reduced only slightly 
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TABLE I: Parameters from analysis of LF-/1SR data in PrOs4Sbi2, Pri_ y Laj / Os4Sbi2 (y = 0.4, 0.6, and 0.8), 
Pr(Osi_ a; Ru :E )4Sbi2 (x = 0.05, 0.1, and 0.6), and PrRu4Sbi2. A: quasistatic K-T relaxation rate obtained from longitudinal- 
field fits to Eq. (2) (see text) except for PrOs4Sbi2. K x : hyperfine enhancement factor (Knight shift) calculated from observed 
low-temperature molar susceptibility (Refs. 29 and 34) and Eq. (4). aw. /i + - 141 Pr Van Vleck relaxation rate (used as measure 
of /i + - 141 Pr coupling strength). r c : correlation time of 141 Pr nuclear spin fluctuations. Experimental values (exp) obtained 
from fits of Eq. (3) to longitudinal field dependence of damping rate A (Fig. 5). Calculated values (calc) obtained from lattice 
sums assuming dipolar coupling (see text). 



Alloy 


A (Ms" 1 ) 


K x 


crvv {fJ-S 1 ) 














exp 


calc 


exp 


calc 


PrOs4Sb 12 a 


0.14" 


19.7 C 


0.26(1) 


0.34 


0.31(2) 


1.31 


Pr .6Lao.40s4Sbi2 


0.21 


16.3 


0.24(1) 


0.22 


0.32(6) 


2.41 


Pro.4Lao.60s 4 Sbi2 


0.21 


14.6 


0.20(1) 


0.16 


0.24(4) 


3.59 


Pro.2Lao.80s4Sbi2 


0.21 


14.7 


0.14(3) 


0.12 


0.2(1) 


5.02 


Pr(Os .95Ru .05)4Sbi2 


0.25 


15.9 d 


0.43(1) 


0.28 


0.19(1) 


1.96 


Pr(Os .9Ruo.i)4Sbi2 


0.24 


12. 2 C 


0.33(3) 


0.22 


0.18(4) 


3.22 


Pr(OS0.4RU0.6)4Sbl2 


0.38 


7.51 c 


0.31(1) 


0.14 


0.31(2) 


7.66 


PrRu 4 Sbi2 


0.22 


6.76 c 


0.18(1) 


0.13 


0.60(6) 


9.19 



°LF-^SR data from Ref. 8. 
''From ZF-^SR data of Ref. 8. 
c Susccptibility data from Ref. 26. 
^Estimated from susceptibility data of Ref. 26. 



for Hl = 10 Oe. Thus the exponential damping is not 
decoupled, which is evidence that A is dynamic rather 
than static in origin. 8 The observed decoupling in Fig. 4 
is associated solely with the behavior of G^ T (A, t), hence 
with quasistatic relaxation by Sb nuclear dipolar fields. 

For any distribution of quasistatic local fields the ZF- 
/xSR asymmetry at long times is expected to approach 
1/3 of its initial value A 23 In Pr-based samples, however, 
the zero-field value of the asymmetry at long times is 
typically much lower than A/3. For PrRu4Sbi 2 this can 
be seen in Fig. 4(c) (circles) (see also Fig. 2 of Ref. 18). 
Like the absence of decoupling, this behavior is evidence 
that the exponential damping is dynamic. 



IV. DISCUSSION 

A. Muon relaxation by fluctuating 141 Pr nuclear 
spins 

Hayano et al. 23 calculated the LF-yuSR and TF-/zSR re- 
laxation rates due to dipole-coupled nuclear spins based 
on the theory of magnetic resonance absorption formu- 
lated by Kubo and Tomita. 36 In their calculation spin 
dynamics arise from muon diffusion in a lattice of qua- 
sistatic nuclear spins, but the treatment applies equally 
to the case where the muon is stationary and it is the nu- 
clear spins that are fluctuating; Lowe and Tse 35 carried 
out an equivalent calculation for relaxation of nuclear 
spins, with essentially the same result. We therefore ap- 
ply the results of Hayano et al. to the present experi- 
ments, after modification to include hyperfine enhance- 
ment of the 141 Pr nuclear spins. This is accomplished by 
replacing the bare 141 Pr gyromagnetic ratio 141 7 by the 
enhanced value 141 7(1 + K), where K is the hyperfine 



enhancement factor (i.e., the 141 Pr Knight shift), 29 ' 34 as- 
sumed isotropic for the tctrahedral 2 Pr site. The damp- 
ing rate A is 

A = g vv ( 3t c 

2 \l + (^H L T c ) 2 



l + [7M- 141 7(l + ^)] 2 ff£r 2 



l + [7 M + 141 7(l + iO] 2 #iT 2 / ' V ; 

where crvv is the high-field muon Van Vleck relaxation 
rate 23 due to 141 Pr dipolar fields (used here as a mea- 
sure of the rms amplitude of these fields), r c is the 
correlation time of the 141 Pr spin fluctuations, and a 
powder average has been taken. The rapid-fluctuation 
limit ervv T c •C 1 is assumed, since otherwise an expo- 
nential damping function would not be expected. 23 This 
relation includes the contributions to A of both longi- 
tudinal Pr spin fluctuations, 35 which are assumed to 
have a common correlation time r c . The field dependence 
of A arises from the fact that the muon spin relaxation 
rate is proportional to the fluctuation noise power at 
the muon Zeeman frequency j^H^; for weakly-coupled 
nuclear spins the dipolar fluctuation spectrum consists 
of broadened peaks centered at 0, 141 7(1 + K)Hl, and 
— [ 141 7(l+if )Hl] , 32 - 35 corresponding to each of the terms 
of Eq. (3). 

From Fig. 5 it can be seen that the field dependence of 
A is generally well fit by Eq. (3). We find, however, that 
the form of this relation does not determine K well com- 
pared to the other fitting parameters. We therefore ob- 
tain K independently from the observed low-temperature 
molar susceptibility x mo i of the Pr ions, 26 ' 27,37 and fix it 
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FIG. 4: (color online). Representative ZF- and LF-/iSR spin 
polarization decay functions, (a) Pro.6Lao.40s4Sbi2. Trian- 
gles: Hl = 5.9 Oe. Diamonds: Hl = 16.1 Oe. Squares: 
H L = 121 Oe. (b) Pr(Oso.9Ruo.i) 4 Sbi2. Triangles: Hl = 
5.3 Oe. Diamonds: Hl = 16.0 Oe. Squares: Hl = 75.6 Oe. 
(c) PrRu 4 Sbi 2 . Circles: Hl = Oe. Triangles: H L = 6.3 Oe. 
Diamonds: Hl = 10 Oe. Squares: Hl = 63 Oe. Curves: 
fits to the damped static K-T function in longitudinal field 
[Eq. (2)] (see text). 



in the fitting. This is done using the relation 29,3 ' 3S 

K = K x = dhfXmol , (4) 

where cthf = 187.7 mole emu -1 is the Pr atomic hyperfme 
coupling constant. The calculated values of K x are given 
in Table I, together with the experimental values of ovv 
and t c obtained from the fits. 



B. Comparison of data to hypernne-enhancement 
model 

Next we compare the experimental values of ovv an d 
t c to those expected from the hyperfine-enhancement sce- 
nario, assuming dipolar couplings between all spins. We 
make the ansatz that the 1 Pr spin fluctuations are due 
to spin-spin interactions within the 1 Pr spin system; in 
analogy with the discussion of muon spin relaxation in 
Sec. Ill A, one would expect significant 141 Pr spin-lattice 
relaxation by 4/ electronic spin fluctuations to result in 
considerable temperature dependence of A contrary to 
experiment. 



We first calculate the "unenhanced" (i.e., using the 
bare 141 Pr gyromagnetic ratio) powder-average /i + - 141 Pr 
Van Vleck relaxation rate er vv from a standard lattice 
sum 23 ' 32 for the candidate (0,^,0.15) muon site 8 (12e 
in Wyckoff notation) in the filled-skutterudite structure 
(space group 7m3). For dipolar 141 Pr- 141 Pr coupling the 
correlation time t® in the absence of hyperfme enhance- 
ment is estimated by the inverse of the unenhanced like- 
spin 141 p r _ 141 p r Van Vleck rate, which is obtained from 
a lattice sum similar to that for cry v . 32 The hyperfmc- 
cnhanced values of these quantities are then 29 ' 34 

<J VY =a%(l + K) (5) 

and 

r c = r c °/(l + K) 2 ; (6) 

in both cases the dependence on K comes solely from the 
hyperfme enhancement of 141 7 (Refs. 29 and 34). 

The calculated values of ovv and r c are given in Ta- 
ble I. For the La-doped alloys ensemble averages of the 
lattice sums have been taken over random 141 Pr site 
locations. 32 The experimental and calculated values of 
<7vv are in rough agreement, but the experimental val- 
ues of t c are systematically smaller than the calculated 
values, sometimes by more than an order of magnitude. 
Correspondingly, the inequality ayyT c <C 1 required for 
our analysis is satisfied by the experimental values of cyv 
and t c but usually not by the calculated values. 

The discrepancy between experimental and calculated 
values of r c is most likely an indication that the as- 
sumption of purely dipolar 141 p r _ 141 p r interactions is not 
valid. Indirect RKKY-like interactions mediated by the 
Pr 3+ intra-ionic exchange interaction 29 may decrease t c 
considerably since they are also hyperfine-enhanced, but 
arc difficult to estimate in PrOs4Sbi2- A brief discus- 
sion of limits on the exchange-mediated 141 Pr- 141 Pr in- 
teraction constant 141 J ex is given in the Appendix, which 
concludes that indirect interactions would account for the 
experimental values of t c (Table I) with only a modest 
Pr 3+ exchange coupling. 

We note that r c is also the 141 Pr NMR signal lifetime 
(spin-echo decay time) T2, so that 141 Pr NMR exper- 
iments would provide an independent measure of this 
quantity. Unfortunately the values of r c from Table I 
(<1 /is) are too short for the NMR signal to be observ- 
able using current spectrometer technology. Moreover, a 
search for the 141 Pr resonance would be difficult because 
of uncertainty in the hyperfine enhancement factor K. 
The mere observation of a 141 Pr NMR signal would there- 
fore significantly modify the conclusions of this paper, a 
fact which might motivate such a search. 

An exchange- mediated 141 Pr-/i + interaction, using a 
Fermi contact interaction between the mediating elec- 
trons and the muon, could also be present. A scalar 
interaction oc I-S^, where I and S M are 141 Pr nuclear 
and muon spin operators, respectively, leads to a muon 
spin relaxation rate of the form of Eq. (3) but with only 
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FIG. 5: Exponential damping rate A vs longitudinal applied field Hl in PrOs4Sbi2 and its alloys, (a): PrOs4Sbi2 (data from 
Ref. 8). (b)-(d): Pri_ y Laj / Os 4 Sbi 2 , y = 0.4, 0.6, and 0.8. (e)-(h): Pr(Osi_ ;c Ru^)4Sbi2, x = 0.05, 0.1, 0.6, and 1. Curves: fits 
to Eq. (3) (Ref. 23). 



the second term in the brackets. 39 Unfortunately the rel- 
ative strengths of dipolar and exchange-mediated inter- 
actions cannot be determined accurately from fits of an 
appropriately generalized version of Eq. (3) to the data. 
Exchange-mediated interactions would, however, be un- 
likely to dominate (hyperfine-enhanced) /i + - 141 Pr dipo- 
lar interactions, since in /-electron metals the magni- 
tudes of dipolar and electron-mediated interactions be- 
tween local electronic moments and muon spins are usu- 
ally comparable. 40 The rough agreement between exper- 
imental and calculated dipolar values of ovv (Table I) is 
consistent with this observation. 

Other sources of uncertainty in the comparison be- 
tween experimental and calculated parameters in Table I 
include (1) the fact that the low-temperature susceptibil- 
ity may not be entirely due to the Van Vleck mechanism, 
leading to error in the calculation of K x , and (2) the 
fact that the calculation of cryv depends on the assumed 
muon stopping site in PrOs4Sbi2. The stopping site has 
not been determined definitively; indeed, high-field TF- 
/iSR 41 suggests that there may be more than one muon 



site. Given these caveats, together with the qualitative 
nature of parameters derived from fitting the /iSR data 
(Sec. IIIB), the agreement between the experimental re- 
sults and the hypothesis of hyperfine-enhanced 141 Pr nu- 
clear magnetism can be regarded as satisfactory. 

Figure 6 gives the Ru and La concentration dependen- 
cies of the zero-field damping rate A(0) just above T c . 
We see that A(0) decreases as the Pr sublattice is diluted 
with La ions, as expected if the dynamic relaxation is 
due to 141 Pr nuclear magnetism. In Pr(Osi„ 2 ,Ru 2; )4Sbi2, 
where the Pr sublattice is not diluted, it can be seen in 
Fig. 6(b) that A(0) generally increases by ^50% as x is 
increased from zero, and decreases again as x — * 1. For 
dipolar coupling A(0) = 5<7 vv t c [cf. Eq. (3)], so that the 
K dependencies of cryv and t c [Eqs. (5) and (6)] can- 
cel and there should be no dependence of A(0) on Ru 
concentration. But this cancellation would not necessar- 
ily hold for an exchange-mediated interaction, which we 
have seen is necessary to understand the short experi- 
mental values of r c . The Ru concentration dependence 
is not well understood, but at least two mechanisms may 
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the muon spin relaxation. 
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FIG. 6: Concentration dependence of zero-field low- 
temperature exponential damping rate A(0) in (a): 
Prj,Lai_yOs 4 Sbi2, (b): Pr(Os a; Rui_ I )4Sbi2. 



be involved: (1) the concentration dependence of AEcef, 
which affects Xmo\ and hence K x [Eq. (4)], and (2) the 
(unknown) concentration dependence of 141 J ex across the 
alloy series. 



V. CONCLUSIONS 

ZF- and LF-^iSR measurements have been carried 
out in the filled skutterudite alloys Pr(Osi_ x Ru x )4Sbi 2 , 
x = 0.05, 0.1, 0.2, 0.6, and 1.0, and (Pri_ y Laj / )Os 4 Sbi 2 , 
y = 0.2, 0.4, 0.6, and 0.8, to clarify the origin of the dy- 
namic damping of the muon spin relaxation observed in 
these alloys. At low temperatures, LF-^SR experiments 
in both alloy series indicate that static local field distri- 
butions and dynamic fluctuations are both involved in 
muon spin relaxation. The temperature and concentra- 
tion dependencies of the muon damping rate A suggest 
that it is due to hyperfine-enhanced Pr nuclear mag- 
netism; the enhancement is responsible for the increased 
141 p r nuc i ear spin-spin interaction strength and conse- 
quent rapid 141 Pr spin fluctuations. Further evidence 
for this picture comes from the field dependence of A, 
which is in reasonable agreement with fits to the model 
of Hayano et al. 23 and calculated coupling strengths (Ta- 
ble I) assuming hyperfme enhancement and exchange- 
mediated spin-spin coupling between 141 Pr nuclei. 

We conclude that hyperfine-enhanced 141 Pr nuclear 
spin fluctuations account for the observed exponen- 
tial damping of the muon spin relaxation function in 
PrOs4Sbi2 and its alloys, and that electronic spin fluctu- 
ations (other than the Pr 3+ electronic response involved 
in hyperfme enhancement) are not directly involved in 



APPENDIX: EXCHANGE-MEDIATED ^Pr- 1 
SPIN INTERACTION 



l Pr 



Consider a 141 p r _ 141 p r S p m interaction Hamiltonian of 
the form 



141 



H 



141 



n_J i 141 n_l 

rtdip + n cx 

141 Jdi P [3(Ii-r 12 )(I 2 -r 12 )/rf 2 



Ii h 



+ 



141 



Jc-J-l ■ I2 



(Al) 



where 141 J d i P = [ 141 7(1 + K)h] 2 /r\ 2 and 141 J cx are the 
interaction constants for 141 p r _ 141 p r dipolar and indirect 
exchange couplings, respectively, and ri 2 is the distance 
between 141 Pr spins Ii and I 2 . A scalar exchange in- 
teraction has been assumed for simplicity. Both 141 Jdi P 
and 141 J cx are hyperfme enhanced. In a simplified model 
(singlet ground and excited CEF-split states), 29,42 the ra- 
tio 141 J7cx/ 141 i7dip is given approximately by 



141 rr /141 7 

dc X j *Jd 



IP 



dip 1 



(A.2) 



where J^. and JJ; are the corresponding interaction 
constants for Pr 3+ interionic exchange and dipolar cou- 
pling. 

To our knowledge has not been determined ac- 
curately in PrOs4Sbi2 or its alloys. Fits to suscepti- 
bility data of CEF models using a molecular-field ap- 
proximation to the exchange coupling 43 yield a molecu- 
lar field constant A = 3.9 mol/emu. This gives J^/hs ~ 
N A nl e X/k B z cS « 0.14 K, where the T = Van Vleck ef- 
fective moment fj, e s « 0.71^b has been assumed and z e s 
is an effective number of nearest neighbors. This estimate 
should probably be considered a rough upper bound. It 
satisfies the criterion 44,45 



(A.3) 



for the absence of exchange-induced magnetic ordering; 
as noted previously, AScef/^b ~ 7 K. There is no 



9 



sign of magnetic ordering in PrOs4St>i2, and many of 
its electronic properties can be accounted for without 
invoking an exchange interaction. The dipolar inter- 
action ^Jjp can, however, be calculated: its maximum 
value is J-L/k B = {gfJ, B ) 2 /k B rl n « 8 x 1(T 3 K, where 



r nn is the nearest-neighbor Pr-Pr distance. A value of 
J7 e x an order of magnitude larger than this would lead 
to 141 J cx / 141 Jdip > 1 from Eq. (A. 2) without violating 
Eq. (A. 3), and could account for the experimental values 
of t c (Table I). 
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